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Fatigue refers to the changes in material properties caused by repeatedly applied loads. It has 
been widely studied for, e.g., construction materials, but much less has been done on soft materials. 

Here, we characterize the fatigue dynamics of a colloidal gel. Fatigue is induced by large amplitude 
oscillatory stress (LAOStress), and the local displacements of the gel are measured through high- 
frequency ultrasonic imaging. We show that fatigue eventually leads to rupture and fluidization. 

We evidence four successive steps associated with these dynamics: (i) the gel first remains solid, 

(ii) it then slides against the walls, (iii) the bulk of the sample becomes heterogeneous and displays 
solid-fluid coexistence, and (iv) it is finally fully fluidized. It is possible to homogeneously scale the 
duration of each step with respect to the stress oscillation amplitude (tq. The data are compatible 
with both exponential and power-law scalings with ao, which hints at two possible interpretations 
in terms of delayed yielding in terms activated processes or of the Basquin law. Surprisingly, we find 
that the model parameters behave nonmonotonically as we change the oscillation frequency and/or 
the gel concentration. 


INTRODUCTION 

Fatigue refers to the changes in material properties 
caused by repeatedly applied loads [1]. Fatigue is quite 
universal and affects a broad range of materials from met¬ 
als [5] to biomaterials such as adhesion clusters of cells 
[3]. In metals, the microscopic events associated with 
fatigue depend on the load and usually weaken the ma¬ 
terial. For example, if loads are above a certain thresh¬ 
old, microscopic cracks form at stress concentrators such 
as the surface, persistent slip bands, and grain bound¬ 
aries [mm. Eventually a crack reaches a critical size, 
propagates, and the structure fractures mi- Here, we 
address the question of fatigue in a soft material namely 
a colloidal gel. Under a constant stress, colloidal gels, like 
many other viscoelastic solids such as foams or emulsions, 
flow only above a critical stress referred to as the yield 
stress ay. This stress-induced fluidization is involved in 
virtually any application of colloidal gels, and it is essen¬ 
tial to correctly measure and control this yield stress to 
ensure the safe and proper use of the material. Yet it has 
been shown in recent years that many phenomena may af¬ 
fect the way a material yields mn]. Slippery boundary 
conditions may induce the solid-body motion of the ma¬ 
terial without bulk fluidization |I2) . Long and complex 
transient regimes may show coexistence of fluid and solid 
phases before a liquid-like stationary state is reached [la¬ 
in particular colloidal gels are prone to so-called “de¬ 
layed yielding” where it may take a very long time for 
a weak gel to yield under a constant stress [T3UTH] . For 
instance, a gel can fail under its own weight after several 
hours or days [nidH]. This “gravitational collapse” has 
been associated to the progressive weakening of the gel 


network due to thermally activated localized rearrange¬ 
ments m or to the interplay between the poroelasticity 
of the gel and gravitational compression [inun]. Delayed 
yielding has also been investigated in gels subjected to a 
constant shear stress imposed by a rheometer, includ¬ 
ing carbon black gels [32 [33j , thermo-reversible protein 
gels [21] or weak gels of polystyrene particles [15]. In 
all cases, the gels are eventually fluidized after complex 
transient regimes. These dynamics have been related to 
activated processes: the applied shear stress cr decreases 
the energy barrier for bond breakage leading to failure 
times that decrease exponentially with a [H] [16]. In 
this work we focus on the interplay between fatigue and 
yielding of carbon black gels under cyclic shear, a phe¬ 
nomenon which is comparatively much less documented 
in soft matter than the case of constant load. 

In a previous paper devoted to large amplitude os¬ 
cillatory shear stress (LAOStress) experiments [37] j 
have first shown that the classical procedure to deter¬ 
mine the yield stress ay displays surprising results. In 
an oscillatory experiment at a constant frequency /, 
where the amplitude Uo of the shear stress is ramped 
up with time t, ay is defined as the stress amplitude at 
which the loss modulus G" overcomes the elastic mod¬ 
ulus G' [38H30]. In carbon black gels, ay depends on 
the velocity of the applied stress ramp [37]: ay decreases 
as the ramp rate increases. This result is inconsistent 
with physical aging, where both G' and ay have been 
shown to increase |3I|. This rate dependence of ay is here 
one of the hallmarks of fatigue. To better character¬ 
ize this fatigue-induced yielding, we performed further 
LAOStress experiments [27], this time keeping the stress 
amplitude constant to Uo = H P^i. An extensive time- 
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resolved analysis of the strain response, of the rheological 
Fourier spectra, and of standard Lissajous plots (stress 
vs strain or shear rate) coupled to simultaneous ultra¬ 
sonic imaging, showed that the gel dynamics involved 
two different timescales, < Tf, where is associated 
with failure at the walls, and Tf with a slower heteroge¬ 
neous fluidization involving solid-fluid coexistence until 
the whole sample becomes liquid-like. The spatial het¬ 
erogeneities observed as the gel slowly fluidizes suggested 
a fragmentation of the initially solid-like gel into macro¬ 
scopic domains eroded by the surrounding fluidized sus¬ 
pension. 

Here, we thoroughly address the influence of the stress 
amplitude and of the frequency, as well as that of the car¬ 
bon black concentration, on such a scenario. In Sect, [j 
we first recall the sample preparation method and the 
specifications of the apparatus for combining standard 
LAOStress rheology and ultrasonic imaging in order to 
probe the local displacement within carbon black gels 
during LAOStress experiments. In Sect, [j we show that 
Tf results from the accumulation of multiple yielding pro¬ 
cesses: the gel first remains solid, at it starts sliding 
against the walls, at u it displays a bulk solid-fluid co¬ 
existence, and at Tf it is finally fully fluidized. While 
the raw times r^,, Tf, and Tf seem to display two different 
behaviours for low vs high stresses (Tq as also observed 
under a constant load in Refs. mMnnna, scaling the 
durations of the successive yielding processes Tu,, Tf, — 
and Tf — Tb with uo allows us to describe the whole stress 
range within a single scaling. These results are further 
discussed in Sect. [] where we compare our data with an 
extension of the delayed failure model m to oscillatory 
shear (Appendix 2) and with the Basquin law. The ex¬ 
tension of the delayed failure model is based on activated 
processes and predicts that the characteristic durations 
decrease exponentially with a^, whereas the Basquin law 
predict a power-law behaviour. Both models fit our data 
well given the small accessible range for uo- Finally, we 
report a surprising behaviour of the parameters of the 
delayed yielding model and the Basquin law when vary¬ 
ing the concentration of the gel and the frequency of the 
oscillations. 

MATERIALS AND METHODS 
Carbon black gel preparation 

Carbon black (CB) particles are colloidal, carbonated 
particles with a typical size range of 85 to 500 nm [33] 
that result from the partial combustion of hydrocarbon 
oils. These particles are widely used in the industry 
for mechanical reinforcement or to enhance the electrical 
conductivity of plastic and rubber materials [33]. When 
dispersed in a mineral oil (density 0.838 g.cm“^, viscosity 
20 mPa.s, Sigma Aldrich), these CB particles (Cabot Vul¬ 


can XC72R of density 1.8 g.cm“^) are weakly attractive 
with interactions of typical strength U ^ SOksT [351 [36] . 
where fcs is the Botzmann constant and T is the abso¬ 
lute temperature. From a dispersed state, the particles 
aggregate to form sample-spanning networks of fractal 
dimension df = 2.2 even at very low concentrations [33| . 
Here, we focus on such gel-forming dispersions at weight 
concentrations of 4, 6 and 8% w/w. 

Standard rheology experiments 

The mechanical properties of colloidal gels are typ¬ 
ically measured using standard rheology experiments. 
Our rheological measurements are performed in a Taylor- 
Couette cell with smooth, polymethyl methacrylate 
(PMMA) walls (height 50 mm, rotating inner cylinder 
radius 48 mm, fixed outer cylinder radius 50 mm, gap 
width 2 mm) with a stress-controlled rheometer (ARG2, 
TA Instruments). The temperature is controlled by a wa¬ 
ter circulation around the Taylor-Couette cell and fixed 
to 25 ± 0.1 °C for all experiments. To ensure a repro¬ 
ducible initial gel state, each measurement is prepared 
using the following sequence of steps: (i) We preshear 
the suspension at 1000 s“^ and then at -1000 s“^ for 20 s 
each to break up any large aggregates, (ii) We let the gel 
rest so that it can restructure by applying a zero shear 
rate for 2 s and then a zero shear stress for 10 s. (iii) We 
probe the viscoelastic properties of the gel by monitoring 
the elastic modulus G' and loss modulus G" for 60 s at 
a very low stress in the linear regime (cto = 0.5 Pa at 
/ = 1 Hz), (iv) The sample is left to rest again for 10 s 
by applying a zero shear stress, which allows for possible 
residual stresses stored during the previous steps to re¬ 
lax [331 [32] ■ Finally, the LAOStress experiment is started 
by applying a sinusoidal stress of amplitude cto and fre¬ 
quency / starting at time t = 0. The reader is referred 
to Refs. [331 HZ] for more details, especially on internal 
stresses and on issues with inertia during LAOStress. 

Ultrasonic imaging under shear 

Several techniques allow for local flow characterization 
under shear, e.g. rheo-optics [351133], particle tracking 
[3nii42| . magnetic resonance imaging [331 HI] or ultra¬ 
sonic velocimetry |45j . Here, we use high-frequency ul¬ 
trasonic imaging |46| . a technique that is fast enough 
to follow the yielding dynamics while being insensitive 
to the opaqueness of the CB gel m- To provide ultra¬ 
sonic contrast the gel is seeded with 1% w/w hollow glass 
spheres (Potters Sphericel, mean diameter 6 /rm, mean 
density 1.1 g.cm“^). This technique is implemented on 
the Taylor-Couette experiment described in Sect. [] A 
linear array of 128 piezoelectric transducers, immersed 
in the water tank and facing the outer cylinder of the 
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Taylor-Couette device, sends short ultrasonic pulses with 
a central frequency of 15 MHz that propagate in a vertical 
plane inclined by about 10° from the normal to the outer 
cylinder. While traveling through the Taylor-Couette 
cell, these pulses get scattered by the hollow glass spheres 
suspended within the CB gel. For each incident plane 
pulse, the backscattered signal, corresponding to the in¬ 
terferences of the various echoes from the hollow glass 
spheres, is recorded by the transducer array, leading to an 
“ultrasonic speckle” signal with 128 measurement lines 
and typically 800 points sampled at 160 MHz. The data 
analysis consists in first processing the speckle signal into 
a beam-formed speckle image S{r, z, t), where r is the ra¬ 
dial direction across the gap (r = 0 being taken at the 
rotating inner cylinder and r = 2 mm at the fixed outer 
cylinder) and z is the vertical direction (z = 0 corre¬ 
sponding to about 10 mm from the top of the Taylor- 
Couette device). Then two successive speckle images are 
cross-correlated in order to get the tangential displace¬ 
ment A as a function of r and z. Here we set the time 
interval between two speckle images to five LAOStress os¬ 
cillation periods. This allows us to map the displacement 
A(r, z,t) of the CB dispersion between successive states 
separated by five oscillation cycles over 32 mm along the 
vertical direction z with a resolution of 250 ^m and over 
the 2-mm gap with a resolution of 100 /im in the radial 
direction r. This temporal and spatial resolution is ideal 
for following the yielding dynamics of the CB dispersion 
under LAOStress without focusing on the intracycle dy¬ 
namics [22]. For full technical details, we refer the reader 
to Refs. PTIHB] . 


Normal force measurements 


In order to measure the normal force N during a 
LAOStress experiment, we switch from a Taylor-Couette 
geometry to a standard cone-plate geometry with cone 
diameter 5 cm and angle 4°. Like in the Taylor-Couette 
geometry, to ensure a reproducible initial state, each 
measurement is carried out by (i) preshearing the sus¬ 
pension at 700 s“^ then at -700 s“^ for 60 s in each di¬ 
rection, (ii) manually shearing the gel until normal forces 
cancel out, (iii) letting the gel rest again for 30 s and 
(iv) measuring G' and G” within the linear regime and 
checking that the normal force remains zero. Finally, 
we start the LAOStress experiment. Calibration experi¬ 
ments have shown that we can neglect the contribution of 
the centrifugal force to the measurements of N provided 
that the shear rate is below 700 s“^. We checked that the 
shear rate never exceeded 100 s“^ during the LAOStress 
experiments. 



/ (s) 



t (s) 


FIG. 1: Time evolution of the elastic modulus G' (blue 
squares) and loss modulus G" (red circles) for a carbon black 
dispersion at c = 6 % w/w under an oscillatory shear stress 
of constant amplitude ctq and frequency / = 1 Hz. (a) For 
a low stress amplitude, ao = 15 Pa, the gel initially dis¬ 
plays a solid-like behaviour {G' > G"), then apparently yields 
{G' = G”) at t ~ 700 s, and thereafter seems to flow like a 
liquid (G' < G”). (b) For a high stress amplitude, ctq = 27 Pa, 
the gel seems to yield and flow like a liquid from the begin¬ 
ning of the experiment. The insets show the corresponding 
time evolutions of the strain amplitude jq. In the soIid-Iike 
regime, the strain remains very low. 70 then increases by two 
orders of magnitudes and finally saturates at some maximum 
value in the fluid-like regime. White squares indicates the 
times at which the displacement field A(r,z,t), as obtained 
from ultrasonic imaging, is illustrated in Fig. 


RESULTS 


Standard rheological measurements 


Figure reports the evolution of the elastic and loss 
moduli of carbon black dispersion at c = 6% under 
LAOStress for both a “low” stress amplitude [cto = 15 Pa, 
Fig 1 'a)] and a “high” stress amplitude [erg = 27 Pa, in 
Fig IT)]. In both cases, when given enough time, the 
sample eventually completely fluidizes: G" > G' while 
both moduli are stationary. However, in the low stress 
case, the gel remains solid until t ^ 700 s whereas, in 
the high stress case, the gel is apparently fluidized from 
the start. To unveil the fluidization processes involved in 
those two different LAOStress experiments, we now turn 
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to ultrasonic imaging. 

Displacement field from ultrasonic imaging 

Figure and movies in the ESI^ show the evolution 
of the tangential displacement held A{r,z,t) recorded 
inside the Taylor-Couette device simultaneously to the 
rheology experiments reported in Fig. with c = 6%. 
One can dehne four successive regimes of material be¬ 
haviour as characterized by different structures of the 
displacement held. First, for 0 < f < t^, A is constant 
and equal to zero across the entire gap. This means that 
from one cycle to the other (or at least every hve cy¬ 
cles since the displacement is only sampled at 0.2 Hz), 
the hollow glass spheres, and thus the gel material, re¬ 
turns into exactly the same position. Together with the 
fact that G' > G", we interpret this as the behaviour 
of a homogeneous solid that adheres to the walls of the 
Taylor-Couette cell. This regime is only observed at low 
stress amplitudes, typically for cto < 20 Pa [see Figs. 0^) 
and[^a) for (Tq = 15 Pa]. For larger erg one directly ob¬ 
serves the following second regime, i.e. one effectively 
has Tuj = 0. For < t < t;,, A is either zero in the bulk 
yet with significant values at the walls [see t = 1450 s in 
Fig. [^a)] or shows a homogeneous non-zero value in the 
bulk [see t = 10 s in Fig.j^b)]. This indicates that the gel 
has yielded near the walls at time t = thus creating 
fluidized lubrication layers at the walls although the bulk 
material remains solid. Consequently, we observe a plug¬ 
like flow with A ~ 0 for low stress but non-zero, constant 
displacement at higher stress amplitudes, which points to 
a global drift of the solid-like gel from one cycle to the 
other and to slippage at the walls. Third, for r;, <t < Tf, 
A is constant over regions spanning a large amount of the 
cell gap yet separated by regions where A shows large, 
fluctuating and apparently random values. As already 
discussed in Refs. [211I17], this is typical of a solid-fluid 
coexistence. Indeed in fluidized zones, the hollow glass 
spheres seeding the dispersion show irreversible displace¬ 
ments from one oscillation cycle to the other due to their 
density mismatch with the CB gel. These displacements 
lead to decorrelation of successive ultrasonic speckle sig¬ 
nals and thus to large erratic variations of A. At low ao, 
this solid-fluid coexistence seems to occur preferentially 
along the vorticity direction [see t = 2780 s in Fig. [^a)], 
whereas at higher erg, we observe a more classical coex¬ 
istence in the radial direction [see t = 230 and 290 s in 
Fig. j^b)]. The fluidized zones progressively grow and 
invade the whole sample until full fluidization is reached 
at t = Tf. Finally, the CB suspension flows like a liquid 
thereafter for t > Tf. 

In Fig. |4j we report the three characteristic times, r^,, 
Tf, and Tf, for different values of the imposed stress am¬ 
plitude of the LAOStress experiment determined using 
the ultrasonic imaging technique. Yielding was only ob¬ 


served in the range ug > 9 Pa. For lower values of the 
stress amplitude, the sample remained solid for at least 
10® s. In Fig. i a vertical, linear path represents the time 
evolution of the gel during a constant stress amplitude 
LAOStress experiment. Then, we can identify a region 
0 < t < Ty,(ao), denoted by (a) in Fig. for which the 
sample is completely solid and adheres to the walls of the 
Taylor-Couette cell. In the region Tw(cro) < t < T{,((Tg), 
denoted by (b), the sample remains solid in the bulk but 
slips along the walls of the cell. We associate this region 
(b) with a fatigue process that will lead to fluidization at 
a later time. In the region Tf,((Tg) < t < r/(crg), denoted 
by (c), the material is heterogeneous exhibiting a coex¬ 
istence of fluid-like and solid-like regions. The solid-like 
zones get progressively eroded by the fluid-like zones un¬ 
til the entire gel is fluidized at t = Tf. A similar erosion 
process has been observed in laponite gels in shear-rate 
controlled experiments [47l |48]. At both low and high 
stresses, r^,, Tf, and Tf seem to follow an exponential law 
of type T ^ situation is complex as one 

can define four different a* related to Tf, and Tf at low 
and high stresses. In Ref. [H], a* at high stresses associ¬ 
ated with Tf is interpreted as the stress barrier necessary 
to be overcome so that the gel yields. In the creep exper¬ 
iments of Ref. m, a similar behaviour is observed and 
a* associated with Tf is modeled in the delayed yield¬ 
ing framework m- In this model, the change of slope 
between the low and high stress regime is attributed to 
a change of the average number of particles involved in 
the weakest link of the gel network mesh. We believe 
that applying the delayed yielding model on Tf is how¬ 
ever inappropriate and that Tf corresponds here to an 
accumulation of multiple yielding processes. 

Normal forces 

On top of this complex spatiotemporal yielding be¬ 
haviour, we observe the presence of normal forces. In 
Fig-i we measure the normal force N in the high stress 
regime during a LAOStress experiment in a cone-plate 
geometry. We observe that G" > G' at t = 0 s: the gel 
directly yields at the wall, which is consistent with the 
high stress experiments in the Taylor-Couette geometry. 
We observe dramatic irreproducible fluctuations of the 
normal force. In our opinion, the value of the normal 
force might be related to residual stress within the gel 
El [37]. However, as soon as the gel partially fluidizes 
in the bulk, for Tf, < t < Tf, we consistently measure a 
negative normal force. This effect is robustly observed 
in all experiments as we varied the frequency / or the 
imposed stress amplitude Cg. 

To understand the effect of the normal forces on the 
erosion process during the solid-fluid coexistence, we as¬ 
sume that deformations within the gel are small, and 
we use the linearized theory of elasticity for the strand 
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FIG. 2: Time evolution of the displacement field A{r,z,t) as measured through ultrasonic imaging during the LAOStress 
experiment shown in Fig. The USV sampling frequency is 0.2 Hz so that A(r, z, t) is probed every five oscillation cycles. 
The rotor is located at r = 0 mm and the stator position is r = 2 mm. A{r,z,t) is coded in linear color levels, (a) For a low 
stress amplitude, ao = 15 Pa, the gel remains completely solid (A ~ 0) for t < Tw — 700 s. For t-u, < t < Tb the gel is fluidized 
at the rotor and the stator but remains solid in the bulk. For Tb < t < Tf the bulk material fluidizes heterogeneously along 
the vertical direction z. For t > Tf the gel is homogeneously fluidized, (b) For a high stress amplitude, no = 27 Pa, the gel 
first displays a homogeneous yet nonzero displacement field indicative of slippage at both walls. The gel subsequently shows a 
heterogeneous fluidization along the radial direction r for Tb < t < Tf and full fluidization for t > Tf. See ESI^ for movies that 
simultaneously display the data recorded by the rheometer and the displacement field A{r,z,t) for both low and high stress 
amplitudes. 



t = 360 s 


t = 620 s 


network which is assumed to be isotropic. The gel is sup¬ 
posed uniform but partially fills the gap. Such a model 
is derived is Appendix 1 and demonstrates how normal 
forces directly affect the stress state within the solid frac¬ 
tion of the solid-fluid coexistence. We find that the pres¬ 
ence of a negative normal force increases the maximum 
principal stress within the gel. Therefore, if the normal 
force is independent of the applied stress ao, the correc¬ 
tion to ao induced by the presence of normal forces only 
shifts Tf to lower stresses but does not affect the scaling. 


Scaling of the characteristic times 

Having excluded the normal forces to explain the dif¬ 
ferent behaviour at low and high stresses, we turn our at¬ 
tention to the duration of each yielding process. We note 
that at low stresses, all the characteristic times have a 
similar value t.^] ^ ^ Tf. However, at high stresses the 

presence of lubrication layer seems to split those char¬ 
acteristic times apart from one another. Given the fact 
that in addition to Tf, we have and t^, we can define 
the duration of each event and compare the low stress 
and high stress regimes, r^, corresponds to the duration 
over which the gel remains entirely solid. Arbw = Th — T^ 
is the duration of the fatigue regime when the gel slips 
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FIG. 3: Time evolution of the normal force during a 
LAOStress experiment for a carbon black gel at c = 6 %w/w 
and / = 1 Hz in a cone-plate geometry at ag = 20 Pa. 
We have repeated the experiment 5 times to check for re¬ 
producibility. (a) Time evolution of the elastic modulus (G', 
solid line) and loss modulus (G", dashed line). The insets 
show the time evolution of the strain amplitude 70 . The be¬ 
haviour of G', G" and 70 is similar to the one observed for 
the Taylor-Couette geometry (Fig.[^b)). (b) Time evolution 
of the normal force N. Inset: average value of N over the 5 
experiments for Tb < t < Tf. The shaded area corresponds to 
the standard deviation. 


at the wall but remains solid in bulk. Aryt, = Tf — Th 
is the duration of solid-fluid coexistence. As shown in 
Fig. i the discrepancy between low and high stresses 
vanishes when those three durations are plotted as func¬ 
tions of the stress amplitude cto- We have therefore fitted 
the data with a single exponential or a single power law 
for each one of the three durations. Given the narrow 
range of experimentally accessible values of ctoj h is diffi¬ 
cult to discriminate between exponential and power-law 
fits (see Table |I] for the values of the best fit parame¬ 
ters in both cases). Still, this scaling allows us to dis¬ 
tinguish three successive yielding phenomena: yielding 
at the wall, bulk yielding and clusters yielding. The gel 
hierarchically yields in time and each successive yield¬ 
ing process requires a higher stress barrier a* or a lower 
exponent a. 



(To (Pa) 


FIG. 4: Dynamic state diagram of the fatigue process of a 
carbon black gel at c = 6 % w/w under LAOStress at a fre¬ 
quency of / = 1 Hz. The characteristic times t = Tw (♦), 
t = Tb (■) and t = Tf (o) are determined from ultrasonic 
imaging and define the boundaries between different states 
encountered during the fatigue process: (a) solid-like state, 
(b) plug-like flow i.e. solid-body displacement with slip at 
the walls, (c) solid-fluid coexistence along the vorticity direc¬ 
tion, and (d) a fluid-like state. Solid lines are exponential fits 
based on Table H] 

DISCUSSION 

Basquin model and delayed failure model 

Gaining insight into the physical principles governing 
our experiments from a theoretical perspective remains 
a challenge. The entire process shows multiple yielding 
time scales, which suggests that many mechanisms are at 
play. The gel first yields at the wall: interaction between 
the wall and the gel should thus be taken into account. 
The gel then yields heterogeneously in the bulk, produc¬ 
ing large clusters that eventually fluidize, which suggests 
that a mean-field theory is inappropriate. Keeping those 
difficulties in mind, and given the fact that no theoreti¬ 
cal approach currently responds to all those criteria, we 
examine two models: the Basquin model and the de¬ 
layed failure model. Both models rely on a mean-field 
approach. We shall therefore apply them on homoge¬ 
neous processes only, such the duration t^ of the initial 
solid regime and that of the fatigue regime A.Tbw 
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FIG. 5: Scaling properties of the characteristic durations as a function of the applied stress amplitude (Jq. (a) Flow state as a 
function of time, (b) Tm, Ari,™ and Ar/;, vs cto- Red solid lines are exponential fits (At = ) to the data whereas blue 

dashed lines in the insets are power-law fits (At = to the same data, (c) Sketch of the gel structure between the stator 

and the rotor. 




























A (10'‘ s) 

a* (Pa) 

B (10' Pa“s) 

a 

'Tw 

9±1 lO’* 

1 .2±0.2 

8±1 10' 

10.5±1.1 

^'^bw 

1 .6±0.2 

5.7±1.1 

l.OiO.l 

3.6±0.5 

Arfb 

1.7±0.2 

6.4±1.0 

3.1±0.3 

3.4±0.4 


TABLE I: Parameters used in Fig. to fit the data by an 
exponential (At = ) and by a power law (At = 

Bor) 

The power-law behaviour of the yielding characteris¬ 
tic times is reminiscent of the empiric Basquin law [32], 
where tj ^ ■ This law can be derived from a fiber 

bundle model [15] . The gel is then represented by a set of 
parallel hbers organized on a regular lattice. When the 
bundle is subjected to an increasing external load, the 
fibers behave like linear elastic springs until they break 
for a given failure load. To recover the Basquin law, Kun 
et al. gsiEnj modified the fiber bundle model and consid¬ 
ered a mean-field approach in which the network strands 
form a random network and fail either due to immediate 
breaking or to aging through damage accumulation. To 
capture damage recovery due to healing of micro-cracks 
or thermally activated rebinding of failed contacts, Kun 
et al. muin] introduced a memory term which scales ex¬ 
ponentially with a characteristic time of the system. In 
this framework, a is directly related to the growth law 
of local damage as a function of the local stress. High 
values of a mean that the material is prompt to accumu¬ 
late damage. Fits of the experiments give = 10.5 and 
Oibw = 3.6 for Tyj and Arf,iu respectively (see Table|I]). The 
value of Ubyj is close to the Basquin exponents found in 
recent creep experiments on carbon black gels {a ^ 2-3) 
P5] . carbopol microgels {a ^ 3-8) [5T| and casein bio¬ 
gels {a ~ 5) [52j but very high when compared to metals 
(a ^ 0.1) [55] and to asphalt {a ~ 0.5) [21||5D]. Soft gels 
are indeed prone to accumulate damage and therefore 
they are much more sensitive to stress than metals. We 
also note that our data show no hint of a critical stress 
as measured in Refs. El [51]. 

To model the exponential behaviour of the characteris¬ 
tic times of the yielding gel, we turn to the delayed failure 
model [16]. This is also a mean field approach, based at 
the micro-scale on Kramers’ transition state theory [55] . 
Here, the particles form percolating strands where n is 
the average number of colloids in the cross-section of the 
strand at its weakest point. When no stress is applied, 
thermal fluctuations are assumed to dissociate or asso¬ 
ciate interparticle bonds at rates kn and Ua respectively, 
with k]j < kA for attractive gels. A constant C with 
units of compliance is also postulated, which converts the 
applied stress a into work Ca with units of ksT. This 
work lowers the interaction potential across each bond. 
Two strand-breaking processes were then proposed: one 
for high-stress case, leading to sequential bond-breaking 
within a strand and an activation stress a* = l/C, and 
one for the low-stress case which allows for re-association 


of bonds (healing) within the strand and an activation 
stress a* = IJnC. Both bond-breaking processes lead 
to an exponential scaling of the time-to-failure. In Ap¬ 
pendix 2, we adapt the theoretical frame developed for 
failure of colloidal gels under a constant load to cyclic 
loading. As for creep experiments, the delayed failure 
model for cyclic loading yields an exponential behaviour 
of the characteristic failure time t™°'^ with respect to oq: 
At = . However, the prefactor A to this expo¬ 

nential is greater than that of the creep, and leads to a 
longer failure time: gels resist better to oscillations than 
to creep for the same activation stress cr* and applied load 
CTo- For Arbw, we find cr* « 5.7. In this model n and C 
are coupled and cannot be measured independently; we 
take values similar to those of Ref. m , ki) = 0.26 s 
kA = 1.5 s“^, l/C = 28.5 Pa and n = 5 as benchmarks 
for further discussion. In the case of r^,, ct* « 1.2 is much 
smaller, implying that the gel is weaker in the wall re¬ 
gion. Keeping C constant, it requires to divide n by ^5 
which means that the CB particle strands develop fewer 
links with the walls as compared to the bulk. Conversely, 
keeping n constant, it requires to divide l/C by about 5 
which means that the CB particles have a much weaker 
interaction with the walls than between themselves. 

Which model is best suited to interpret the present 
experiments remains an open question. The fiber bun¬ 
dle model can solely support longitudinal deformation 
which allows for studying only loading of the bundle par¬ 
allel to fibers so that its use to model our experiment 
performed under shear is at best qualitative. Moreover, 
the Basquin law stipulates that the gel will always fail 
even at very low stresses. Experimentally, it seems that 
below roughly 9 Pa the gels remains solid, even though 
it cannot be excluded that it fails after some very long 
time exceeding experimental time scales. In the delayed 
yielding model, a divergence of the time-to-failure can 
be explained through the dominance of healing in the 
competition between healing and rupture of strands. A 
predictive capability for this critical stress is not yet avail¬ 
able. We also note that u* = 6±1 Pa roughly matches 
the stress below which it seems impossible for fatigue to 
eventually fluidize the gel. 

Influence of the frequency of the LAOStress 

So far we have discussed the measurements and the in¬ 
terpretation of the scaling behaviour of the characteristic 
times with (Tq . We now focus on two parameters that are 
likely to modify the scaling parameters of Tf with cto, 
namely the CB concentration c and the oscillation fre¬ 
quency /. In general, / and c do not qualitatively affect 
the fatigue scenario described in Fig. [^ Moreover, in the 
high stress regime, the Tf{ao) data are well fitted both 
by an exponential and by a power law. Figure [^ shows 
the values of a* and a as a function of / and c. 
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As for exponential fits, the stress barrier a* increases 
sharply with c. This seems reasonable since a concen¬ 
tration increase reinforces the backbone structure of the 
gel. On the contrary, the exponent a found from power- 
law fits does not display any systematic evolution with c. 
Although we probe higher <jq for higher c, the suscepti¬ 
bility of the gel appears to be independent of c within the 
fiber-bundle model interpretation of the Basquin law. 

The major advantage of oscillatory experiments com¬ 
pared to creep experiments [22l [23] is that LAOStress 
is frequency-resolved. Figure |^a) shows that a* first 
increases with / at low frequencies, then reaches a max¬ 
imum around /=0.4±0.15 Hz and finally decreases for 
higher frequencies. The delayed yielding model does not 
predict any frequency dependence of a* (see Fig. 0^) 
in Appendix 2). Yet, we know that the gel rebuilds it¬ 
self on time scales that are below 1 s [52]. Therefore 
decreasing the frequency should leave more time for the 
gel to heal between two successive oscillation cycles and 
make the gel more resistant to stress. This could explain 
the decreasing curve cr*(/) at high frequencies. On the 
other hand, a displays a minimum around /=0.4±0.15 
Hz while the Basquin model also does not predict any 
frequency dependence since it only depends on the num¬ 
ber of cycles. Yet the increase of a at high frequencies is 
qualitatively consistent with the previous healing argu¬ 
ment that leads to increased susceptibility to stress as / 
increases. We note that Z. Shao et al. ]SS] have also ob¬ 
served a puzzling, nonmonotonic frequency dependence 
in the yielding process of flocculated microgel suspen¬ 
sions. Currently, we do not have any sensible argument 
to interpret the nonmonotonic variations of a or a* with 
/ as reported in Fig. [^ and more specifically their low- 
frequency behaviour. 

CONCLUSION 

We have investigated the fatigue scenario in LAOStress 
experiments on carbon black gels. Yielding was only ob¬ 
served within experimental time scales when the applied 
stress amplitude exceeded a critical value (cto > 9 Pa for 
c = 6% w/w). Above this applied stress, we have de¬ 
termined three characteristic times for delayed yielding 
under LAOStress based on ultrasonic imaging: within 
the time interval 0 < t < the gel remains solid, within 
Tw < t < Tf, the gel yields at the walls but remains solid 
in the bulk, and within Tb < t < Tf , we observe solid- 
fluid coexistence associated with negative normal forces. 
Beyond Tf, the sample flows like a liquid. Those char¬ 
acteristic times scale differently with the applied stress 
(To at low and high stresses. However this difference van¬ 
ishes when we scale the characteristic durations of each 
process r^,, Tb — Tf — Tb with (Tq. Such a scaling al¬ 
lows us to distinguish three successive yielding processes: 
yielding at the wall, bulk yielding and clusters yielding. 



FIG. 6: Evolution of (a) the stress barrier a* and (b) the 
power-law exponent a with the oscillatory stress frequency / 
for a carbon black gel at c = 4 (■), 6 (•) and 8% w/w (♦). 
Lines are included to guide the eye. 

Those durations are functions of the applied stress and 
can be empirically fitted by an exponential, although a 
power law cannot be ruled out due to the limited range 
of the stress amplitude. To gain insight into those multi¬ 
ple yielding processes, we have used the delayed yielding 
model extended to LAOStress (t ~ exp(— (To/ct*)) and 
the Basquin model (t ^ o'q°‘). Both models fit reason¬ 
ably well the data. This analysis shows that each succes¬ 
sive yielding process requires a higher stress barrier, a* 
or a lower exponent a. According to the Basquin model, 
our gels are prompt to accumulate damage and are much 
more sensitive to stress than asphalt [33] or metals [53] . 
According to the delayed failure model, the carbon black 
particles have a weaker interaction with the walls than 
between themselves which justifies that the gel yields first 
at the wall and then in bulk. 

The yielding process also depends on the concentration 
of the gel and on the applied frequency in a complex man¬ 
ner that remains to be fully understood. Another inter¬ 
esting line of enquiry, left for future study, is the effect of 
the particle-wall interactions on the yielding at the wall. 
Presumably, weak or repulsive particle wall interactions 
would reduce the time-to-failure at the wall, while 
strong attractive particle-wall interactions could prolong 
Tu,, or even re-localize initial yielding to the bulk. In 
the same spirit the influence of surface roughness, which 
has already been addressed in carbon black dispersions 
submitted to a steady shear stress [53], should also be in¬ 
vestigated and modelled in the case of LAOStress. 
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Finally, there is some similarity between the cyclic fa¬ 
tigue of the colloidal gel and the very high cycle fatigue 
(VHCF) [56] observed for metal alloys. When subjected 
to a cyclic stress far below the yield stress, microcracks 
develop dispersed within the crystal grains. These mi¬ 
crocracks multiply and grow until they reach the grain 
boundaries. This triggers an avalanching crack growth 
and ultimate macroscopic failure. Similarly to carbon 
black gels, the time-to-failure for VHCF in metals decays 
exponentially with the stress amplitude. Hence, despite 
the very different microstructures of metals and colloidal 
gels, there are striking similarities in their fatigue be¬ 
haviour that presumably derive from thermally-activated 
microscale processes governing irreversible damage. 


then gives 


(1 - iy)crii - U<7-22 = 0 
(733 - i/(crii -I- 0 - 22 ) = 0 


(la) 

(lb) 


where v is the Poisson modulus of the gel. The stress 
tensor becomes 


0 

a{t) —N{t) 0 

0 0 


( 2 ) 


The principal tensile stress is obtained by computing the 
maximum eigenvalue of a, which gives 


Appendix 1: Normal forces model 


= - 


2{l-u) 


N{t)- 


ll- 412(1 -v) 
4(1 - ^)2 


N‘^{t) -I- fT^(t). 


To potentially understand the origin of the normal 
forces during the solid-fluid coexistence, we assume that 
deformations within the gel are small and we use the lin¬ 
earized theory of elasticity for its strand network. The 
strand network is further assumed to be isotropic. We 
take the gel to be in a state of uniform, plane strain, with 
zero deformation in the z direction of the Taylor-Couette 
geometry or the normal direction of the cone-plane geom¬ 
etry. We also neglect the curvature of these geometries. 
Thus, we introduce a rectangular coordinate frame, with 
base vectors { 61 , 62 , 63 }, so that e} is the shear direction 
and 6*2 is the shear gradient direction. From the plain 
strain geometry, we have cri3 = (T31 = (T 23 = o ’32 = 0. 
The remaining stress components depend on the bound¬ 
ary conditions. When the gel region fills the gap of the 
geometry and adheres to the walls, the boundary condi¬ 
tion is 

(712 = a{t), 

while we have en = 622 = £33 = 0 inside the gel, where 
6 denotes the Euler-Almansi strain tensor. The above 
yields the well-known simple shear solution 

0'12 = 0-21 = cr(t), 

while all other stress components vanish for infinitesimal 
strains. The principal tensile stress then becomes identi¬ 
cal to the applied shear stress crp{t) = 

During solid-fluid coexistence, the gel only fills a frac¬ 
tion of the gap, and the boundary conditions on the gel 
become 

ai2 = a{t), a22 = -N(t), 

where N(t) is the normal force observed in experiments. 
From the geometry, en = £33 = 0 inside the gel. 
Hookei’i^s law for plane strain and the equilibrium equa¬ 
tion for an isotropic, linear elastic material in plain strain. 


Therefore, the presence of a negative normal force in¬ 
creases the principal stress on the gel and vice versa. 

Appendix 2: Delayed failure model for cyclic loading 

The delayed failure model is a mean field approach 
that describes the dynamics of the fraction of remaining 
strands, x, in a yield surface of a colloidal gel subject to a 
creep experiment at constant stress a(t) = (Jq. At start, 
the gel is completely solid and x =1. The dynamics is 
then driven by type ordinary differential equation (ODE) 
of birth/death type [TS] : 

dx 

— =-KD(x,a{t))x + Ka{1 - x); (3) 

Kr){x,aQ) = Kq exp( "^‘^° ) is the dissociation rate of 
the strands and Ka is their association rate. depends 
on local parameters such as kA and ki), respectively, the 
association and dissociation rate of the individual col¬ 
loids, n the average number of colloids in the cross sec¬ 
tion of a strand at its weakest point, the temperature T, 
the range of the attraction between colloids 6, and the 
initial area density of strands pq. Indeed, the exponential 
behaviour oi Ku is related to an activated process where 
the stress barrier necessary to break a strand composed 
of n colloids in its cross section is 1/nC = poksT/S. The 
complete expression of the amplitude ol Kjj can be found 
in Ref. [TO] . 

Instead of looking at a creep experiment where cr{t) = 
CTq we now consider LAOStress experiments where the lo¬ 
cal stress on the strands is taken to be a(t) = ctoI sin 27r/t| 
in Eq. (|^. Taking sin27r/t instead of |sin27r/t| would 
lead to Kjj ~ 0 when the sin function becomes negative, 
which does not correspond to the actual situation where 
the sample experiences the same dissociation rate inde¬ 
pendently of the sign of the applied shear stress. We then 
discretize x into x„i = x(m/f), the fraction of the strands 
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in the gel at the cycle, m being an integer. Inte¬ 
grating Eq. ^ over the m*'' loading cycle with respect 
to time, while assuming that x changes only marginally 
during each cycle, gives a recurrence relation 


Hm+i)// fnCao \ 

Xm+i-Xm = -Xm / E'q cxp -|sin 27 r/t| dt 

Jm/f \ / 

+ ^(1 - a^m) 

= ~X^T (+ ^(1 - Xm). ( 4 ) 

/ \ Xm J J 


with initial value a;o = 1, and where 


. 1/2 

T(a) = 2 / exp(a| sin27rs|)ds. (5) 

do 

Using again that Xm+i—Xm is small, the recurrence equa¬ 
tion @ can be approximated by a differential equation, 
now viewing m = /t as a real number: 


dx 

dt 


= -A'.T 

dm 


f nCaQ 


' + Ka{ 1 - x). 


( 6 ) 




This ODE describes how distributed damage evolves 
under cyclic loading conditions. To the knowledge of the 
authors, no numerical solution exists to Eq. (§. For 
a sufficiently large Ka > 0, that is if ruptured strands 
are assumed to re-associate, the numerical solutions show 
that X converges toward a non-zero equilibrium value re¬ 
lated to the ratio Ka/Kq- no failure can be predicted 
with certainty. Still, experimentally we always observe 
a binary behaviour of the yielding process: either erg is 
too small and the suspension remains solid, x = 1 over 
some seriously long time, or the suspension eventually 
completely fluidizes, x = 0. We therefore restrict the fol¬ 
lowing discussion to the fluidization process and assume 
that Ka = 0. 

The numerical solution of the ODE in Eq. ([^ with 
Ka = 0 gives the strand density remaining attached x 
as a function of various parameters such that Kq, the 
oscillation frequency / and the amplitude of the applied 
stress <7q. As shown in Fig.[^ x decreases from 1 to 0 as a 
function of time: the gel gets completely fluidized. In this 
model, provided that the oscillation inverse frequency is 
small compared to the yielding time, the frequency does 
not influence the yielding process [see Fig. [^a)]. erg sets 
the timescale of the fluidization process [see Fig. [^b)]. 
We extract two characteristic times from this numerical 
resolution: the time Tfiuid = t(x = 0) at which x = 0 
and the characteristic yielding time Tqsc defined as the 
inverse slope of the tangent to x at t=0. Tosc is the time 
that matches the closest the definition of the experimen¬ 
tal duration r^,, the time necessary for the gel to yield 
at the wall or Axbw, the duration necessary for the gel 
to yield in bulk. Indeed, in both those regimes and in 


FIG. 7: Resolution of the delayed failure model for 
LAOStress. (a) Density of strands x remaining solid under 
LAOStress as a function of time for different frequencies at a 
constant applied stress amplitude no = 6 Pa. (b) Density of 
strands x remaining solid under LAOStress as a function of 
time for different applied stresses at / = 1 Hz. (c) Character¬ 
istic times obtained by solving Eq. ([^ numerically 

(□), and Eq. ([^ analytically (—). (-) is the 

analytical solution of the delayed failure mode l for the corre¬ 
sponding creep experiments given in Eq. (101. The parame¬ 
ters of the model are displayed in Table [II 


n C (Pa-^) 

1/nC (Pa) 

ko (s kA (s 

Axbw 5 0.033 

6 

0.26 1.5 


TABLE II: Parameters used in the numerical resolution of 
the delayed failure model for LAOStress and chosen so as 
to roughly match the values of Ati,w in Fig. [^ Note that 
different sets of parameters can give the same solutions to the 
ODE. Indeed n and C are coupled and as long as the stress 
barrier IjnC is constant, the solutions are identical. 


the region of interest x is close to I. Note that Tfiuid 
should not be compared to r/ as the model assumes an 
homogeneous yielding of the strands; experimentally we 
know that this is not the case. The model is thus only 
optimal to describe the experiments close to x = I when 
the gel is solid and Fig. [^c) only shows Tq^c vs ug. 

Although no analytic solution is available for Eq. §, 
we can find a solution for the limit x close to 1. The 
solutions of Eq. (§ crucially depends on the properties 
of T. For strong gels m, like carbon black gels, failure is 
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only observed at experimental time-scales if nC(T q ^ 1. 
Therefore, the asymptotic behaviour of T for large (Tq 
is pivotal. For large values of ag, the integral defining T 
is dominated by the short interval near s = 1/4 where 
the sine function has its maximum. Consequently, we 
may use a Taylor expansion for sin 2tts around s = 1/4 
and change the limits of the integral to ±oo, this way 
obtaining 
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T(a) « 2 


rjA‘ 


1 - 


(27rs — 1/4)^ 


ds = \ —e“, 
na 

(7) 


for a > 1. Using this approximation with Eq. (§ gives 


( 8 ) 

Given that the solution is reminiscent of an exponen¬ 
tially decaying function, we obtain a characteristic time- 
scale for fatigue failure of strong gels with Ka = 0, simply 
by taking 


= nCWo^ nCao > 1 . 

. 

This can be compared to the corresponding time in a 
creep experiment, Tcreep, which was previously found to 
be [T51 

Tcreep = ^ (nC'cro)“^e“"'^‘^° . (10) 

The numerical solutions to Eq. Q with Ka = 0 compare 
excellently with the analytical solution of Eq. Q when 
sampled at the frequency /, nCao 3>1 and x close to 1 
[see Fig.[7](c)]. 

To conclude, the delayed yielding model shows that the 
characteristic failure time varies exponentially with (Tq for 
both LAOStress and creep experiments: Tosc Kreep 
exp(—CTo/(T*). The intimate relation between creep and 
cyclic fatigue originates from the common microscopic 
process governing failure at low stresses in the case of 
this relatively strong colloidal gel. cr* = 1 /nC represents 
the stress barrier necessary to break the gel. In the case 
of oscillations, a* is frequency-independent. 
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